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We evaluate the quasinormal modes of massless Dirac perturbation in a Schwarzschild black hole 
surrounded by the free static spherically symmetric quintessence by using the third-order WKB 
approximation. The result shows that due to the presence of quintessence, the massless field damps 
more slowly. The real part of the quasinormal modes increases and the the absolute value of the 
imaginary part increases when the state parameter Wq increases. In other words, the massless Dirac 
field decays more rapidly for the larger Wq. And the peak value of potential barrier gets higher as 
... |fe| increases and the location of peak moves along the right for fixed Wq. 

, 

O \ PACS numbers: 04.30.Nk, 04.70.Bw, 97.60.Lf 

^ " ^ . Keywords: Quasinormal modes; Dirac perturbation; WKB approximation. 

(N 



o 

(N 



> 



X 
S3 



I. INTRODUCTION 



Quasinormal modes of black holes play a important role in the black hole theory. They are regards as the "charac- 
teristic sound" of black holes. The concept of QNMs was firstly pointed out by Vishveshwara[H in calculations of the 
scattering of gravitational waves by a black hole. In recent decades, many people have done a lot of work of the quasi- 



I ' normal modes of black holes [1|-[T3]. They have investigate the perturbation such as: gravitational perturbationplj . 



^jPj scalar perturbation electromagnetic perturbation [13|i Dirac perturbation|14|- |17j, a nd so on. There are n iany 
methods to calculate the QNMs of the pert urbation, for example, WKB approximationp^-(20j . the "potential fit" [2ll |. 
and the method of continued fractions[2^, etc. In this paper, we use the third WKB approximation. 
There are three major motivations to study the QNMs in detail. First, QNMs may provide valuable help to identify 
black hole parameters 2^ 23|. Second, an important aspect of QNMs studies has been related to the AdS/CFT conjec- 
ture. By computing QN frequencies in AdS spcectime, we can obtain a prediction for the thermalization timescale in 
the strongly coupled CFT. It is argued that string theory in AdS space is equivalent to CFT in one less dimension [2^. 



And the third is the relation between QNMs and black hole area quantization[26|-[28 
Recently, V. V. Kiselev^] has considered Einstein's field equations for a black hole surrounded by the static spher- 
ically symmetric quintessential matter and obtained a new solution that depends on the state parameter Wq of the 
[ — quintessence. And Chen and Jing evaluated the quasinormal frequencies of massless scalar field perturbation fsoj 
, around the black hole which is surrounded by quintessence. In other papers, we have investigate the gravitational 
L| ■ perturbation 31 1 and the electromagnetic perturbation ^32*]. In this paper, we consider the massless Dirac field in this 



situation. 

II. DIRAC EQUATION IN THE SPACE TIME OF SCHWARZSCHILD BLACK HOLE SURROUNDED 

BY QUINTESSENCE 

The metric for the spacetime of Schwarzschild black hole surrounded by the static spherically-symmetric 
quintessence is given bv[30l| 

ds"^ = ^fdf + f-^dr^ + r^{de^ + sin^edcf)^), (1) 

with 
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where M is the black hole mass, Wq is the quintessential state parameter, c is the normalization factor related to 
Pq = ~§ ^3(1+^,) ; and pq is the density of quitenssence. 

The massless Dirac equation in the black hole spacetime can be written asfssll 

[7''e;^(a^ + r^)]* = o, (3) 

where is the inverse of the tetrad e° defined by the metric g^i, , 

gtiy = 'riabeti"'eJ' (4) 

with rjab ~ diag(— 1, 1, 1, 1) being the Minkowski metric. 

7° is the Dirac matrix, and is the spin connection which is given by 

^, = \b\l'Vae,,;, (5) 

where ehu-^fj, = df^ebu — ^'^u^ba is the covariant derivative of Cbu, and F"^ is the Christoffel symbols, 
we take the tetrad as 

e"^ =diag{^/f,-^,r,r sin 0). (6) 



F^ can be expressed as 



Fo = i/'7'7°,ri = 0,r2 = iy77S',r3 - i (sin 0/77173 + cos 07V) (7) 



Then, the Dirac equation ([3]) becomes 

^/f dt \dr r Af dr J r 09 2 r sin 9 dtp 



7° 5* ^ 1 / a 1 1 7^ 9 1 , 73 a^- 



Define 

5- = /-^$, (9) 

Eq. ((HI) becomes 



-y// at \or r J r oO 2 rsmO dtp 



Define the tortoise coordinate as 



Introducing the ansatz 



dr 
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with spinor angular harmonics 



/,+ ^ I V 2/+1 



l+^+m Y^n— 1/2 1 ^ 

lH-rn'+l/2 )' ^' = ^+2^' 

21 + 1 



(10) 



(11) 



*=fST'^^"'^''''M^-'"'' (12) 



( / '+a-"' ym-l/2 \ 
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equations [10] can be written in the form 

The cases for (+) and (—) in the functions and can be put together giving 

f£ + iuj'-V^)F = 0, (14) 



with 



drl 



^''""WTZ+Jf -/). {fork = i + \, „„<i3 = i + i), (16) 



V. - ^(WV7-^| + /), [for + (17) 

The potentials Vi and V2 are supersymmetric partners derived from the same superpotential. We shall concentrate 
just on Eq. [14] with potential Vi in evaluating the quasinormal mode frequencies for the massless Dirac field by the 
third WKB approximation in the following sections. 

III. MASSLESS DIRAC QUASINORMAL FREQUENCIES 

In this section we evaluate the QN frequencies for the massless Dirac field in the space time of the Schwarzschild 
black hole surrounded by quintessence using the third WKB approximation. WKB method was firstly developed by 
B. F. Schutz and C. M. Will fisl at the second order, later, S. Iyer and C. M. Willp^ developed the method to the 
third order and R. A. Konoplya[20j extended it to the sixth order. The accuracy of the WKB formula is better with a 
larger multipole number I and a smaller overtone n. The formula for the complex quasinormal frequencies a;(for the 
third order WKB method) is 



a;2 



where 





(77+ 188a2) 



lOOa^ 



(19 + 28q^) 



Vr 



(4)- 



2304 \ Vo" 



(67 + 68a2) 




(19) 



and 



r,=r,(rp) 



(20) 



Take M — 1, c = 0.001 and M — l,c — Q for our calculation. And c = means there is no quintessence. Using the 
third-order WKB approximation, we can get the solutions as the table 1 and table 2 show(for ^ n < k), where 
I is the angular harmonic index, n is the overtone number, oj is the complex quasinormal frequencies, Wq is the 
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FIG. 2; Variation of the effective potential for the massless Dirac field in the Schwarzshild black hole surrounded by quintessence 
with r for |fc| = 5, c = 0.001 and 3t«, + 1 = -0.01, -1.7, -1.99. 



quintessential state parameter. 

The effective potential V{r, fc), which depends on the absolute value of k and Wq, is in the form of a barrier. In Fig. 
1, it shows the variation of effective potential with when Wq = — ^ and c — 0.001. The variation of the effective 
potential with r which is respective to the quintessential state Wq parameter for fixed |fc| = 5 and c — 0.001 is shown 
in Fig. 2. 

TABIE I:The quasinormal frequencies of massless Dirac perturbations in the Schwarzshild black hole[l3] without 
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quintessence (c=0) . 
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TABIE II: The quasinormal frequencies of massless Dirac perturbations in the black hole surrounded by quintessence 
for |fc| = 2, \k\ = 3, |fc| = 4, = 5 and c = 0.001. 
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FIG. 3: The relationship between the real and imaginary parts of quasinormal frequencies for the rnassless Dirac perturbations 
in the background of the black hole surrounded by quintessence for fixed c = 0.001 and no quintessence for c = 0. x and *, 
respectively, refering to the value with quintessence and without quintessence. 



IV. DISCUSSION AND CONCLUSION 



The quasinormal modes of a black hole present complex frequencies, the real part of which represents the actual 
frequencies of the oscillation and the imaginary part represents the damping. 

From Fig.l we see that the peak value of potential barrier gets higher as |A;| increases and the location of peak moves 
along the right for fixed Wq. From Fig. 2, it shows that when the absolute Wq increases, the peak value of potential 
barrier gets lower. Using the third WKB approximation, we calculate the lowly decaying modes frequencies for 
k = 1,2, 3, 4, 5. The complex quasinormal mode frequencies for the massless Dirac field are listed in Table I(without 
quintessence) and Table II(exist quintessence, and c=0.001). We plot the relationship between the real and imaginary 
parts of quasinormal frequencies with the variation of u>q(for fixed c = 0.001), compared with the situation without 
quintessence. From Fig. 3 we can find that for fixed c(unequal to 0) and / the absolute values of the real and imaginary 
parts decrease as the absolute value of the quintessence state parameter Wg increases. It means that when the absolute 
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value of Wq is bigger, the oscillations damp more slowly. And the absolute values of the real and imaginary parts of 
quasinormal modes with quintessence are smaller compared with those with no quintessence for given I and n. That 
is to say, due to the presence of quintessence, the oscillations of the massless Dirac fields damp more slowly. 
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